Single-order-parameter description of glass-forming liquids: A one-frequency test 
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Thermo- viscoelastic linear-response functions are calculated from the master equation describing 
viscous liquid inherent dynamics. From the imaginary parts of the frequency-dependent isobaric spe- 
! >— ■ ^ ■ cific heat, isothermal compressibility, and isobaric thermal expansion coefficient, we define a "linear 

dynamic Prigogine-Defay ratio" At p {uj) with the property that if At p (lo) — 1 at one frequency, 
O , then Atp(w) is unity at all frequencies. This happens if and only if there is a single-order-parameter 

description of the thermo-viscoelastic linear responses via an order parameter (which may be non- 
P_| ' exponential in time). Generalizations to other cases of thermodynamic control parameters than 

, temperature and pressure are treated in an Appendix. 

<n : 

i 1 I. INTRODUCTION 

A famous result of classical glass science is that if a glass-forming liquid is described by a single order parameter, 
the Prigogine-Defay ratio is unity. Recall that, if Ac p is the difference between liquid and glass isobaric specific 
heat per unit volume at the glass transition temperature T g , Akt the similar liquid- glass difference of isothermal 
compressibilities, and Aa p the difference of isobaric thermal expansion coefficients, the Prigogine-Defay ratio II is 
definedi^ by 

a 

g: n= Ac ] Anr_ 

T g (Aa p ) 2 

\l ■ Both experimentally and theoretically the following inequality was established early on: II > 1,1,2,3,4,5,6,7,^9 The yast 
majority of reported IPs are significantly larger than unityi^, and for decades the consensus has been that - possibly 
with a few exception o 11 ^ 12 - a single order parameter is not sufficient for describing glass-forming liquids. 

The above picture came about as follows. The glass transition is a freezing of configurational degrees of freedom. If 
kinetic aspects are ignored, the glass transition has the appearance of a second order phase transition in the Ehrenfest 
7-H ' sense with continuity of volume and entropy, but discontinuity of their thermodynamic derivatives. If the "order 
\& I parameters" (numbers characterizing the structure) are denoted by Zi,...,Zk, these are frozen in the glass phase. 
In the equilibrium liquid phase the order parameters depend on pressure and temperature, a dependence that is 
determined by minimizing Gibbs free energy G(T,p, zi, ...Zt)- dG/dzi = 0.— In the classical papers by Davies and 
Jones from the 1950's^ it was shown that within this framework one always has n > 1 and that II = 1 if there 
is just a single order parameter. The simplest dynamics for the approach to equilibrium for the order parameters, 
z.i = —AidG/dzi, imply exponential decays for each order parameter for small perturbations from equilibrium. Under 
this assumption, in the case of just one order parameter any quantity relaxes exponentially to equilibrium after a small 
disturbance. Exponential relaxations of the macroscopic variables are seldom observed, although there is evidence 
that relaxation is often intrinsically exponential so that the non-exponential behavior comes from a distribution 
of exponential relaxationsj 13 ' 14 d 5 This reflects dynamic heterogeneity, but historically the fact that macroscopic 
relaxations are non-exponential was often taken as confirming the conventional wisdom that one order parameter is 
rarely not enough. - Besides the reported Prigogine-Defay ratios being almost always significantly larger than unity, 
a further classical argument for the necessity of more than one order parameter is the following'^ If structure were 
characterized by a single order parameter, glasses with the same index of diffraction would also have all other physical 
properties identical, which is not the case. This does not rule out the possibility that a single order parameter 
is sufficient to describe the linear thermo-viscoelastic response of the viscous liquid phase, however, and it is this 
possibility we shall inquire into in the present paper. 

Because kinetic aspects cannot be ignored, the glass transition is not a phase transition. Thus there is no exact T g 
and strictly speaking the n of Eq. |T]) is not well defined. This is because the properties of a glass depends to some 
extent on the cooling history; moreover, extrapolations from the glass phase are somewhat ambiguous because glass 
properties change slightly with time. Thus the changes in specific heat, etc., found by extrapolating glass and liquid 
specific heats from below and above, respectively, to TgH are not rigorously well defined. This weakens the generally 
accepted conclusion that in the vast majority of cases there must be more than one order parameter. 
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The problem of making the Prigogine-Defay ratio rigorously well defined is solved by referring exclusively to linear- 
response experiments on the equilibrium viscous liquid phasei 16 ' 18 ' 19 In viscous liquids thermodynamic coefficients are 
generally frequency dependent, c p = c p (lu) etc. The high-frequency limits reflect glassy behavior where relaxations 
do not have enough time to take place, corresponding to "frozen" structure. Thus for the equilibrium viscous liquid 
phase at any given temperature T one redefines Eq. (fT]) to be the "linear Prigogine-Defay ratio" give n 16 i 18 ' 19 by 



^ _ [c p (uo — > 0) - c p (uo -> oo)] [kt(u — y 0) - kt(uj -> oo)] 
T[a p (u> — y 0) — a p (to — > oo)] 2 

Descriptions of viscous liquids and the glass transition by means of one or more order parameter were actively 
discussed in the 1970'si 16 ' 18 ' 19 ' 20 ' 21 ' 22 Since the II of Eq. (JTJ) is not well defined and since there are no measurements 
of the linear Prigogine-Defay ratio of Eq. ([2]), the question of one or more order parameters remains open . 23 ' 2 In our 
view, there are a number of reasons to take seriously the possibility that some glass-forming liquids may be described 
by a single order parameter: 1) In early computer simulations Clarke reported Prigogine-Defay ratios close to one for 
a 216 particle Lennard- Jones model of Argon^i This is not a highly viscous liquid, but recent extensive simulations 
by Sciortino and co-workers of viscous liquids confirmed this scenario by showing that structural relaxations for 
small temperature changes are basically controlled by a single parameter (e.g., the volume) i 26 ' 27 ' 28 ' 29 2) Experiments 
monitoring loss peak frequency and loss maximum of the Johari-Goldstein dielectric beta relaxation process upon 
temperature jumps revealed a striking correlation between these two quantities, a correlation which is difficult to 
explain unless structure is parameterized by a single order parameter^ In these experiments the relaxation of beta 
process properties is controlled by the alpha relaxation time. This is to be expected because the alpha time is the 
structural relaxation time, but it is highly non-trivial that alpha and beta relaxations correlate in equilibrium viscous 
liquids, as shown recently by Bohmer and co-workersj^I this may be taken as a further indication that a single 
order parameter controls the behavior. 3) Experiments studying dielectric relaxation under varying temperature and 
pressure conditions showed that the shape of the alpha loss peak as quantified by the stretching exponent f3 depends 
only on the loss peak frequency^ 2 - This could be a consequence of there being just one order parameter, determining 
both alpha relaxation time and stretching. 4) Finally, we would like to mention theoretical works suggesting that 
the Prigogine-Defay ratio may be unity in some casesi 23 ' 33 ' 34 - More speculative, in our opinion it would not be too 
surprising if there is just one order parameter for liquids where time-temperature superposition applies accurately, 
because in these liquids linear relaxations appear to be particularly simple with a generic w -1 / 2 high-frequency decay 
of, e.g., the dielectric and mechanical alpha loss peaks i 35 ' 36 

The wide-frequency measurements required to determine II of Eq. (J^J) will be very difficult to perform in the 
foreseeable future. There are, in fact, no measurements of all three required frequency-dependent thermo-viscoelastic 
linear-response functions. Work at our laboratory indicates that such measurements are possible^ but only over a 
limited frequency range. This motivates a search for an alternative to the linear Prigogine-Defay ratio. In this paper 
we introduce a "linear dynamic Prigogine-Defay ratio" (where the double prime indicates the imaginary part of the 
corresponding linear-response function): 



A Tp (uj) 



ToKH) 2 



(3) 



Based on the description of viscous liquid dynamics in terms of Markovian inherent dynamics in Sees. II and III 
we prove that At p (lo) — 1 if and only if the linear Prigogine-Defay ratio (Eq. l[2"])) is unity. This happens if and 
only if there is a single, generally non-exponential order parameter. Moreover, it is proved that if At p (ui) = 1 at one 
frequency, this is true at all frequencies. In section IV we investigate the matter in a thermodynamic approach and 
prove that in this framework the existence of a single order parameter implies that At p (u>) — 1 at all frequencies. In 
the Appendix it is shown that At p (uj) is just one out of a family of 4 linear dynamic Prigogine-Defay ratios larger 
than or equal to one, which equal unity if and only if there is a single (generally non-exponential) order parameter. 



II. RESPONSE MATRIX FOR MARKOVIAN INHERENT DYNAMICS 



Adopting the energy landscape approach to viscous liquid dynamics^ 8 - we model the liquid (henceforth: "system") by 
the inherent dynamics consisting of jumps between the potential energy minima l 29 ' 39 ' 40 ' 41 This description is believed 
to be realistic in the highly viscous phase where there is a clear separation between the vibrational time scales (in 
the picosecond range) and the (alpha) relaxation time scale - the time scale of interest here. If the system has N 
potential energy minima, an ensemble of similar systems is described by a vector of probabilities P = (Pi, Pi...., Pjv) 



3 



where P n denotes the fraction of systems vibrating around energy minimum n. If G„(T, p) denotes the vibrational 
Gibbs free energy of minimum n, the P-dependent Gibbs free energy is defined^ by 

G(T,p,P) = Pn(G n (T,p) + k B T In P„) . (4) 

n 

Via the standard thermodynamic identities the P-dependent entropy S and volume V are given by 



S(T,p,P) = -Y,P n (^ + k B \nP n \ 



V(T, P ,P) = (B) 

n 

It is convenient to introduce the notation X = (T, — p) for the controlled ("input") variables and Q = (S, V) for the 
induced responses ("output"). Thus 



Qa = —qx~> a=l,2. (6) 

The equilibrium probability distribution is found by minimizing G(T,p, P) under the constraint P n — 1; this 
leads^ to the well-known canonical probabilities (where G(T,p) is the Gibbs free energy of Eq. ^ evaluated at the 
equilibrium probabilities): 



Pn( T ,P) = ex P ■ ( 7 ) 



k B T 

Inherent dynamics consist of transitions between minimal The inherent dynamics are modelled as a Markov 



process; thus the probabilities obey a standard master equatio: 



,43 



Pn = ^W nm P m . (8) 

m 

The rate matrix W and the equilibrium distribution P£ q both depend on T and p, and for all n one has 



J2w nm (T,p)P^(T, P ) = 0. (9) 

m 

In the following the variables T and p are assumed to vary slightly from their average values To and po as arbitrary 
externally-controlled functions of time, resulting in small changes of S and V. Examining perturbations around the 
reference state Xo = (To, — po) it is convenient to introduce the notation W° m = W„, m (To,po) and P° = P T j 9 (To,po). 
If 5 denotes perturbations from the reference state, P n (t) = P° + 8P n (t) etc., Eq. ([5]) implies to first order 



SPn = £ (KmSPm + 5W nm P^) . (10) 
m 

Since 5(£ m W nm P%) = (Eq. ©), one has £ m <W nm P° = - J2 m W° m 5P^. When this is substituted into Eq. 
(fTTJjl we get 



SP n = Y,Km(6Pm-6P%) . (11) 

m 

Next, we expand 5P^ in terms of 5X a . Equations {5| and imply 44 (where (3 = 1, 2 and f^i is the Kronecker delta 
symbol) 
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dX k B T \dP, 

Thus to lowest order 



ami- _ i f°Q?_ Qf + SflkB , (12) 



8P eq = 



p: 







kBTo 



dQl --Qi + k B )5X 1+ (^-Q 2 )5X 2 



dP m ^ a 11 \ dP, 



(13) 



Inserting this into Eq. (fTTj) and utilizing Eq. ((9|) we arrive at the following "equation of motion" for the probabilities 
when temperature and pressure vary slightly as arbitrary functions of time (SXp = 5Xp(t)): 

SP n = J2< m [sPm rVE^I^^ I 1 

Consider now a harmonic linear perturbation of the system. In relaxing systems like glass-forming liquids the 
ordinary thermodynamic response functions are frequency dependent and complex. Writing T(t) = To + Kc(5T e lu)t ) 
and p(t) = po + Re(5p e luJt ), the oscillations of S and V are similarly described by S(t) = So + Re(SS e lu)t ) and 
V(t) — Vo + Ke(SV e luJt ). If SP n and SXp denote the (complex) frequency-dependent amplitudes, in steady state Eq. 
(fj"4"f implies that 

£ « - iorfta) SP n = ^ J2 WLP^SXp . (15) 

n m,/3 m 

For > we solve Eq. (|15p by defining (where / is the identity matrix) 

leading to the following equation for the frequency-dependent complex probability amplitudes 



6P n = Y^Anm^^SXp. (17) 
m.p 

In order to calculate the thermo-viscoelastic linear-response functions we first define = fx^- This 2x2 matrix 
gives the instantaneous (glassy) linear response because it determines the variations of S and V when the probabilities 
P n are frozen^ Next, we expand the frequency-dependent complex amplitudes SQ a to get 

Here and henceforth, for any function / the notation f(To,po) for a variable that also depends on the probabilities P 
signifies that / is evaluated at the equilibrium probabilities (Eq. (0). When Eq. (jTTJ) is inserted into Eq. (TTg|) , we 
get 



where 



8Q a = Y,J a p(u)SX l3 , (19) 




J af3 (u) = ^+E^(7 1 o,Po)4 m H§ £ (To,Po). (20) 
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Since A nm (uj — > c>o) = it follows that J a p(ui — > oo) = J^g, justifying the notation. 

We proceed to show that the 2 x 2-matrix J a p{oj) is symmetric. Introducing the matrix Y nm defined by 



F nro = (P n °)-^ n ° m (P°)^ (21) 

the detailed balance requirement implies that Y is symmetric ! 43 ! 45 For later use we note that Y is negative 
semidefinitei^ 3 - For all vectors x one has (x\Y\x) < 0. Moreover, equality applies if and only if x oc (P )?; the 
latter property of Y expresses "ergodicity," i.e., the assumption about the master equation that all states are con- 
nected by some sequence of transitions. We define a diagonal matrix R by 



Rum = [Pn) k 5 nm (22) 

and note that W° = RYR^ 1 . This implies that 



Ma) = -r\r(W° -iiuiy^R 2 = —J— (R(Y - iujI^R- 1 ) 1 RYR = —^—-R(Y — iwI)~ l Y R . (23) 
k B T k B T k B T 

Since Y and R are both symmetric, it follows that A(oj) is symmetric. Thus by Eq. (f2TH) the 2x2 matrix J a p{u)) is 
symmetric. 

The J a p(uj) matrix determines the frequency-dependent thermo-viscoelastic linear response defined as follows: If 
Cp(uS) denotes the isobaric heat capacity per volume, a p (oj) the isobaric thermal-expansion coefficient, and kt(uj) the 
isothermal bulk compressibility, the complex frequency-dependent coefficients <5T, 8p, 5S, and SV are related by the 
following matrix (the symmetry of which reflects the symmetry of J a p{ui))\ 



( 5S\ _ y ( Cp(uj)/T a p (u) 
[SV J - Vo { a p {u) k t (lo) 

Whenever there is time-reversal invariance, the (dc, i.e. zero-frequency) symmetry implied by a Maxwell relation 
translates into (ac) Onsager reciprocity, a well-known resultj 43 i 45 




III. LINEAR DYNAMIC PRIGOGINE-DEFAY RATIO 



We proceed to calculate the linear Prigogine-Defay ratio of Eq. @ from the expression for the frequency-dependent 
linear response functions. Since j4 nm (oo) = 0, we find by inserting Eq. (|20[) into Eq. ^ (with T = Tq and where 
dS/dP n = dS/dP n (T ,po) and dV/dP n = dV/dP n (T a ,p a ) but u (T ,po)" is left out for brevity) 



(0) (0) 

n = 2 • (25) 

(0) 

Next we reason in a way analogous to that of Davies and Jones^ based on the Cauchy-Schwartz inequality. This 
inequality is the well-known mathematical theorem that if a real symmetric matrix B is positive or negative semidef- 
inite, the following applies: For any vectors x and y one has (x|B|?/) 2 < (x|B|a;)(y|B|y), and for non-zero x and y 
equality applies if and only if a number 7 exists such that x — jy £ N(B) where AT(B) is the kernel (null space) of 
B. The matrix A(0) is positive semidefinite and the kernel of A(0) is the one-dimensional linear subspace spanned 
by the vector (1, 1)4^ This implies that II > 1 and that II = 1 if and only if there are constants 7 and c such that 
the following equations apply for all n (compare Refs.—' 5 -'^): 



l^ (To ' Po) = 7 H; (To ' Po) + c - (26) 

For reasons given in the next section this situation is referred to as the single-order-parameter case. 

As mentioned in the introduction, the linear Prigogine-Defay ratio of Eq. ^ is difficult to measure. Instead we 
propose to consider the "linear dynamic Prigogine-Defay ratio" defined in Eq. ([3]). When Eq. ([2H)l is inserted into 
this we find 
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A Tp (u) = 




Because the matrix A" (to) is negative semidefinite for all u> > 0^ the Cauchy- Schwartz inequality implies that 
At p (uj) > 1. Again, equality applies if and only if a number 7 exists such that the vector with n'th component 
dV/dP n — 'ydS/dPn is in the kernel of the matrix A"(ui), which is the one-dimensional space spanned by the vector 
(1, 1)4£ Thus for all u> > the equation At p (lo) = 1 is mathematically equivalent to Eq. (f2T))) . and At p {lo) = 1 is 
equivalent to II = 1. In particular, if At p {lo) is unity at one frequency, At p {uj) is unity for all frequencies - and this 
happens precisely when 11=1. 



IV. WHY II = 1 CORRESPONDS TO A SINGLE ORDER PARAMETER 



We define the order parameter e by 



Combining Eq. (|28| with Eqs. (fl8| and (|26[) we find that in any linear experiment (thus calculating to lowest order) 
S and V are functions of T, p and 8e (J^ — J21): 



SS(t) = 6e(t) + Jfi5T(t) - JfS6p(t) 

5V(t) = 7 Se(t) + J£5T(t) - J%8p(t) . (29) 

The situation described by these equations is more general than the single-order parameter model of Prigogine and 
Meixner,— because Eq. (|29j) allows for an order parameter with non-exponential dynamics. The common physics, 
however, are that besides T and p one single number determines both entropy and volume. In the present approach 
e might well have a memory of the thermal prehistory. 

By the definition we shall adopt here, the single-order-parameter situation applies whenever Eq. (I29p is obeyed for 
some variable Se. In last section we proved that if At p (u>) = 1 applies at one frequency, then Eq. (f29|) follows. We 
proceed to show that conversely, if Eq. (|29p applies for some order parameter s, then At p (lo) = 1 for all lo. To prove 
this, consider a situation with periodically varying temperature and pressure fields. In steady state SS, SV, and Se 
vary periodically with complex amplitudes 8S(uj), 6V(u>), and 6e(u>). According to Eq. ([29|) these amplitudes are 
given by (for any u>) 

8S(uj) = Se(uj) + J^ST(uj) - J™6p(uj) 

5V(u) = jde(uj) + J£8T(uj)-J™5p(uj). (30) 
In the case where only temperature varies, comparing to Eq. (|24j) shows that the following two equations apply: 



V (6S(u)\ (5e{lo) 

T C ^ = {S7W)) P = {ST& 



V a p (u) = (^4) = 7 fM) (3D 



5T(u;)J p '\ST(oj) /p 



Since the J°°'s are real numbers, we have 



= a (uj). 



(32) 



Similarly, if only pressure varies Eqs. (|30|) and (I24p imply 
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V a p (u>) 



8S{lo) 
5V{uj) 



Sejuj) 
Sp(uj) 

UpM 



+ j. 



22 i 



(33) 



thus 



7 <( W ) = (34) 

Eliminating 7 by combining Eqs. (f32|) and (|34|) yields At p (w) = 1 for any uj. Note that in the case of one order 
parameter the imaginary ( "loss" ) parts of all three response functions are proportional. 

The above line of reasoning can be repeated for different choices of input and output variables. Altogether there 
are 4 different natural (linear) dynamic Prigogine-Defay ratios as demonstrated in the Appendix; here it is also shown 
that the requirement of positive dissipation implies that no dynamic Prigogine-Defay ratio can be smaller than unity. 



V. CONCLUDING REMARKS 



The original Prigogine-Defay ratio of Eq. (JTJ) is not rigorously well defined. The "linear" Prigogine-Defay ratio 
of Eq. {2} is well defined, but requires measurements of thermo-viscoelastic response functions over many decades 
of frequency. There are still no methods for measuring a complete set of these response functions. Hopefully such 
measurements are possible in the near future, but realistically they will initially only cover a few decades.— This 
motivates the one-frequency criterion proposed and developed mathematically in this paper: If the linear dynamic 
Prigogine-Defay ratio is unity at one frequency, this quantity is unity at all frequencies and a single-order parameter 
description applies. Conversely, if a single-order parameter description applies, the dynamic Prigogine-Defay ratio 
is unity at all frequencies. If this happens, the imaginary parts of the three linear response functions of Eq. (|24j) 
are proportional. As shown in the Appendix, these results all generalize to the 3 other natural choices of input and 
output variables. - These are strict mathematical statements; in practice one cannot determine whether or not the 
dynamic Prigogine-Defay ratio is precisely one. What one can do is test how close to unity this quantity is. We expect 
that a single-parameter description is a good approximation whenever the dynamic Prigogine-Defay ratio is close to 
unity. Note that if there are several relaxation processes with well separated time scales, it is possible that some of 
these are well described by a single order parameter, though not perfectly, whereas others are not. In this case the 
frequency-dependent one-parameter test developed here in principle would be useful although, as mentioned already, 
it will probably not be possible in the foreseeable future to determine all three required response functions over wide 
frequency ranges. 

In the present paper we referred to the At p = 1 situation as that of a single, generally non-exponential order 
parameter. As emphasized by Goldstein,^ however, due to a mathematical equivalence it is really a matter of taste 
whether one prefers instead to refer to a situation of several order parameters obeying a mathematical constraint 
(in our case Eq. ([26)) ). It is not possible a priori to estimate how restrictive or unlikely it is that Eq. (|26[) applies; 
only experiment can settle this question. It should be noted, though, that there is a simple physical interpretation 
of the single-order parameter case: In the approximation where each inherent state is regarded as a potential energy 
minimum with a harmonic potential, via Eqs. and (JT]) Eq. (|26p is equivalent to V n — a\E n + «2 where V n 
is the inherent state volume and E n the inherent state energy (potential energy minimum). Thus there is a single 
order parameter if and only if inherent state volume correlates perfectly and linearly with inherent state energy. It is 
expected that this condition is obeyed to a good approximation if and only if the dynamic Prigogine-Defay ratio is 
close to unity. 
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APPENDIX: GENERALIZATIONS TO OTHER CONTROL VARIABLES 



Standard thermodynamics give rise to a number of (dc) linear-response coefficients. If the thermodynamic variables 
of interest are T,p, S, V, there are 24 coefficients of the form (da/db) c with a,b, and c chosen among T,p, S, 
These coefficients form 12 pairs that are trivially related by inversion [e.g., (da/db) c = l/(db/da) c ]. As is well known, 
the 12 coefficients are not all independent, but related by various Maxwell relations (summarized below). There 
are the following 8 basic linear-response coefficients (where the specific heats are per unit volume and the last three 
coefficients have no generally accepted names): 



Isochoric specific heat: cy = 

Isobaric specific heat: c p = p^ff^ 

Isothermal compressibility: ht = — y(jj^j 

Adiabatic compressibility: us = — y^j^j 

Isobaric expansion coefficient: a p = y{^^j = — "^(ff) 
"Adiabatic contraction coefficient:" as = — p 

"Isochoric pressure coefficient:" (3v = (fr) v = (w). 

"Adiabatic pressure coefficient:" (is = (jr) = (^) 



v 

T 

S 

P- 

s . 



(A.l) 



Consider harmonically varying scalar thermal and mechanical perturbations of equilibrium for a small volume 
element. "Small" here means that its linear dimensions are much smaller than both the heat diffusion length and 
the sound wavelength corresponding to the frequency under consideration, implying that the perturbations may 
be regarded as homogeneous over the small volume element. Let ST(lu), 5p{oj), Ss(uj), Sv(lo) denote the complex 
amplitudes of perturbations varying with time as cx exp(io^) where we define the intensive variables v = V/Vq and 
s = S/Vq (Vo is the equilibrium volume). If small perturbations around temperature To are considered, following 
Eq. (jA.l [) one defines the following complex frequency-dependent linear-response quantities (where according to the 
so-called correspondence principle^ all thermodynamic relations survive and the Maxwell relations become Onsager 
reciprocity relations^): 

Isochoric specific heat: cy(w) = T ^ g^fy ^ 



Isobaric specific heat: c p (lo) = To y ^ 

Isothermal compressibility: kt{uo) = ~ 0jf$p^ J 



( Sv(u) 



5s(u 
8p(u 



(A.2) 



T 
V 



Adiabatic compressibility: ks(co) 

Isobaric expansion coefficient: a p (cu) = ^ st(u) ) 

"Adiabatic contraction coefficient:" ag(cj) = — ( ^["] ) 

"Isochoric pressure coefficient:" (ivi^) = {w0^j = ( Jv0) ) 

"Adiabatic pressure coefficient:" Ps{u) = g?pfy ) g = ( 8v(u) ) 

The general procedure now works as follows. Let (5X(ui), SY (u)) denote the amplitudes of two periodically varying 
thermodynamic variables considered as control variables and (S Z (u) , 6W (u)) the remaining two periodically varying 
variables. The relationship between the two sets of variables generally takes the form 

( 5Z(lo) \ ( aiiH oia(w) \ / 5X(lo) \ 

{ SW(lo) J ~ \ a 21 {u) a 22 {io) ) \ SY(cu) J ' [A - 6) 



In the case of a single order parameter e - compare to Eq. ([50]) - we can write (where the 2 x 3-matrix is real and 
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frequency independent) : 



V sw(u>) j \ a % af 2 b 2 j y 



(A.4) 



Comparing to Eq. I|A.3[) . when Y does not vary we get: 



f =>• » ; x = 7- = 7xy , (A.5) 
J aula;) c>i 



and when X does not vary: 



oiaH = aj§ + h (p$) x | a ;' 2 (co) = b_ 



,- = 7xy ■ (A. 6) 



In principle there are 6 different choices of control variables. Below we treat the 4 natural cases where the input 
variables consist of one from the (5, T) "energy bond" 51 ' 52 and one from the (— p, V) energy bond, and similarly for 
the output variables. In each case the signs are chosen to make the response matrix symmetric. Via Eq. (IA.5j) and Eq. 
(|A.6p in each of the 4 cases the 3 imaginary (loss) parts of the relevant linear response functions are proportional. The 
4 cases are detailed below, where the explicit amplitude frequency dependence is left out for simplicity of notation. 
1. Control variables ST and — Sp . 

This case was considered in the main paper, but is included for completeness: 



Ss \ = / c p (uj)/To a p (uj) \ / ST 
Sv J \ a p (uj) h t {uj) ) \ -Sp 



(A.7) 



Applying Eq. (|A.5[) and Eq. (|A.6[) to this we get 



?oapM _ _ _ 4( w ) 



and thus 



c' p ' (w) p dp (ui) ' 



A Tp (cj) ee fi Lp> = 1. (A.9) 



2. Control variables Ss and Sv : 



To KH]' 



ST \ = ( T /c v (uj) -l/a g (fij) \{5s 
-Sp J [-l/a s (w) 1/k s (uj) J\Sv 



(A.10) 



Applying Eqs. (1A.5|) and (1A.6|) to this we get 



(To/cyH) (l/«5(w)) 

and thus 



— 7SV = -- — ; -—77 (A. 11) 



AsyH (ro/cvH) (1/^))' = L (A. 12) 

V/as^))" 1 



3. Control variables Ss and — <5p : 
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ST 
—5v 



T /c p (lo) -1/0 s {w) \ f Ss 



Applying Eq. (|A.5|) and Eq. (|A.6[) to this we get 



(r /cp(w)) 



= 7s P = 



and thus 



Asp(^j = — = V - 

(V&M)" 



(A.13) 



(A. 14) 



(A.15) 



4. Control variables <5T and Sv 



Ss \ _ ( cv(u>)/T -0v(v) \ ( 6T 
-Sp J I -/3y(w) -1/ktM J V 



Applying Eq. (|A.5|) and Eq. (|A.6|) to this we get 



(A.16) 



(A.17) 



and thus 

Atv 



_ < vH(i/« r H)" = 
To KM] 2 



For each case we showed above that, if there is a single order parameter, the dynamic Prigogine-Defay ratio is unity. 
Below we proceed to prove that, conversely, if one of the dynamic Prigogine-Defay ratios is unity at some frequency 
then this dynamic Prigogine-Defay ratio is unity at all frequencies and there is a single order parameter. In particular, 
if one dynamic Prigogine-Defay ratio is unity at one frequency, all dynamic Prigogine-Defay ratios are unity at all 
frequencies. 

In the periodic situation the dissipation is proportional to lm[ST(u))Ss* (u>) — Sp(uj)Sv* (o;)]i 48 ' 53 The requirement of 
positive dissipation implies that all 4 dynamic Prigogine-Defay ratios are larger than or equal to unity; this follows 
by considering the special case of in-phase input variables in which case it is easy to show that the determinant of the 
imaginary response matrix must be non-negative to have positive dissipation. If one of the dynamic Prigogine-Defay 
ratios is unity at some frequency u, the determinant of the imaginary part of the corresponding response matrix 
is zero. This implies that there is an eigenvector of the imaginary response matrix with zero eigenvalue. Thus the 
dynamic Prigogine-Defay ratio is unity if and only if a thermodynamic cycle exists with zero dissipation. It follows 
that, if in one of cases 2-4 the dynamic Prigogine-Defay ratio is unity at frequency u>, then At p (uj) = 1. This implies 
that At p (uj) = 1 at all frequencies (Sec. Ill), and that there is a single order parameter e such that S and V are 
given by Eq. (|29p. These equations are easily rewritten to give the required output variables in terms of the input 
variables and e. This implies that the relevant dynamic Prigogine-Defay ratio is unity at all frequencies, and that all 
other dynamic Prigogine-Defay ratio are also unity at all frequencies. 
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